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I . Introduction 

The  problem  of  wave  propagation  along  a horizontal  wire  parallel  to 
the  earth's  surface  was  first  investigated  nearly  half  a century  ago  by 
Carson,  Pollaczek,  and  Haberland  [1-5]  in  a low  frequency  approximation 
suitable  for  application  to  overhead  power  lines.  Since  that  time,  various 
refinements  to  the  low  frequency  theory  have  been  made  [6-8]  of  which 
excellent  summaries  can  be  found  in  the  books  [9,10],  However,  it  was 
eventually  desired  to  determine  higher  frequency  characteristics  of  such 
wires,  in  order  to  accurately  account  for  the  earth's  influence  on  elevated 
antennas,  as  well  as  to  describe  transient  processes  on  overhead  power  lines. 
Thus,  about  twenty  years  ago,  more  rigorous  analyses  of  the  problem  began 
to  appear  [10-25].  (See  also  the  additional  bibliography  at  the  end  of  the 
reference  list.) 

With  few  exceptions  [10,  11,  14,  21-24],  the  excitation  of  the  wire 
has  not  been  considered,  and  even  these  are  limited  to  a delta-function 
voltage  generator  [10,  14,  22-24],  plane-wave  incidence  [21],  or  a purely 
formal  treatment  of  a magnetic  ring  current  [11].  The  remainder  of  the 
work  has  sought  to  find  modal  solutions  for  the  source- free  case,  by 
postulating  a current  distribution  I exp{iwt  - Tz} , where  T is  an  unknown 
propagation  constant  along  the  wire,  evaluating  the  resultant  primary  fields 
and  secondary  (induced  by  the  earth)  fields,  and  enforcing  the  boundary 
conditions  at  the  wire  to  obtain  an  eigenvalue  equation  for  T . In  its 
most  rigorous  formulation  [10-17,  20],  this  equation  is  an  implicit  tran- 
scendental equation  accurately  solvable  only  by  numerical  methods.  Only 
under  suitable  approximations  [12-15,  17-20]  can  a closed  form  expression 


for  F be  obtained. 


The  modal  approach  in  essence  is  an  interpretation  of  the  electro- 
magnetic characteristics  of  the  wire  from  a waveguide  point  of  view.  This 
structure,  being  both  open  and  lossy,  can  possess  a quite  complicated  set 
of  modal  properties  [26 ] . We  may  expect  the  discrete  modes  to  be  finite 
in  number,  and  thus  require  an  additional  set  of  continuous  or  radiation 
modes  to  form  a complete  modal  set  in  terms  of  which  an  arbitrary  field 
of  a source  may  be  expanded.  Indeed,  we  may  find  that  improper  leaky  modes 
may  exist  on  this  structure,  as  has  been  claimed  by  some  researchers  [16,25]. 
Alternatively, it  is  possible  that  the  mode  set  of  this  structure  does  not  form  a 
complete  set.  This  is  a well-known  occurrence  in  the  theory  of  nonselfadjoint 
operators  [27].  In  fact,  one  could  call  into  question  the  physical  reality 
of  the  modes  themselves  on  the  basis  of  the  fact  that,  if  T = a + i3,  with 
a>0,  the  assumed  current  distribution  on  the  wire  must  become  infinite  at 


Z = -co. 


In  view  of  such  difficulties,  the  most  rigorous  approach  to 
the  problem  would  seem  to  be  the  analysis  of  a wire  of  finite  length  and 
with  finite  excitation  by  an  integral  equation  approach  typical  of  antenna 
theory.  That  this  approach  possesses  its  own  difficulties  is  shown  by  the 
work  of  [10,23]  in  which  a number  of  approximations  is  made  until  the 
current  distribution  along  the  wire  becomes  essentially  that  of  a standing 
wave  with  a wave  number  identical  to  that  of  the  discrete  mode  found  by 
the  aforementioned  infinite  wire  analyses  [15,17].  Apparently,  then,  the 
modal  viewpoint  is  not  unrelated  to  the  antenna  theory  viewpoint  for  an 


Even  on  a lossless  structure  for  which  a=  0,  one  could  argue  that  the 
assumed  current  distribution  violates  the  outgoing  radiation  condition  at 
z = -<»,  unless  some  source  were  postulated  to  exist  there. 
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antenna  of  reasonable  length.  On  the  other  hand,  one  may  argue  from  the 
analysis  of  an  imperfectly  conducting  coaxial  line,  for  example,  that, 
given  an  arbitrarily  long  line,  the  current  on  the  wire  would  eventually 
be  dominated  by  radiation  currents  instead  of  "modal"  currents  which 
usually  decay  in  an  exponential  fashion  [28].  Thus,  a discrete  modal 
representation  is  useful  only  if  the  length  of  the  wire  is  such  that  the 
current  along  the  entire  length  of  the  wire  can  be  adequately  described 
by  discrete  "modal  currents"  for  a prescribed  excitation  of  finite  extent. 

The  foregoing  discussion  points  up  the  need  for  a rigorous,  general 
analysis  of  the  modal  representation  of  the  fields  on  infinite  thin-wire 
structures  parallel  to  the  air-earth  interface.  Thus,  we  seek  a general 
expression  for  the  excitation  of  the  modes  of  propagation  along  the  wire 
by  an  aperture  field  of  arbitrary  form  in  the  transverse  plane.  Some 
further  numerical  investigation  of  the  behavior  of  the  discrete  modal 
propagation  constants  is  given  in  Section  III.  We  also  present  some  field 
plots  for  the  discrete  modes  for  a variety  of  wire  parameters  are  presented, 
and  the  differences  between  the  modes  are  examined  in  some  detail.  A 
knowledge  of  the  modal  pecularities  is  vital  to  the  ability  to  excite  one 
of  them  with  any  efficiency. 

II.  Excitation  of  a horizontal  wire  above  the  ground 
by  a given  current  distribution 

Let  us  consider  the  problem  of  exciting  an  infinite  thin  wire 
parallel  to  and  above  the  earth  by  a prescribed  set  or  electric  and 
magnetic  current  sources  3®xt  and  J®xt  in  the  region  outside  the  wire 
as  shown  in  Fig.  1(a).  The  height  of  the  center  of  the  wire  above  the 
earth  is  h,  the  radius  of  the  wire  is  a,  the  region  above  the  earth 


"rr: ; '-■? 


_• — : 


r: 


6 

( z > 0)  is  taken  to  be  free  space  with  permittivity  eQ  and  permeability 

UQ,  while  the  earth  (z  < 0)  is  considered  to  have  a complex  dielectric 
2 

constant  e = n where  n is  a complex  refractive  index.  A time 

dependence  of  exp(-iu)t)  will  be  assumed. 

The  solution  of  this  problem  can  be  applied  to  some  other  related 

problems  as  well.  If  the  wire  is  semi-infinite  (z  > 0)  as  shown  in  Fig.  1(b), 

the  effect  of  the  sources  can  be  represented  approximately  in  terms  of  the 

aperture  field  E°  produced  by  the  sources  at  z = 0.  This  would  be  an 

idealization  of  the  physical  situation  sketched  in  Fig.  1(c),  where  the 

wire  is  driven  by  a voltage  generator  of  some  kind.  Following  Harrington 

[29],  we  may  replace  the  sources  and  aperture  field  by  a perfect  (electric) 

conductor  at  z = 0 with  an  adjacent  <etic  current  sheet  E°  x a 

J t z 

located  at  z = 0+.  This  in  turn  is  equivalent  to  removing  the  conductor 
and  imaging  the  magnetic  current  sheet  (as  well  as  the  wire  and  the  medium) 
to  give  a total  current  J°6(z)  = 2E°xaz6(z)  in  the  aperture  plane 
exciting  an  infinite  wire  over  the  earth  (Fig.  1(d)). 

The  wire  will  be  assumed  to  be  thin  compared  to  a wavelength  and  to 
the  height  above  ground: 

A = k^a  « 1;  A « H = k^h 

where  k^  = w/iT^eT  is  the  wavenumber  of  free  space.  Under  these  conditions, 
the  equivalent  total  current  I(z)  may  be  taken  to  be  the  average  over 
the  circumference  of  the  wire  of  the  H field  component  multiplied 
by  2tra  [14].  The  average  value  of  the  Ez  field  component  at  the 
circumference  of  the  wire  will  then  consist  of  three  parts.  The  first 
contribution  is  from  the  induced  (unknown)  equivalent  current  I(z) 
[11,14,15,24]: 


7 


|X 

1 

1 

© 

(o) 

APERTURE  PLANE 

EXTERNAL 

SOURCES 

1 

1 

1 

€0  ' Mo 

1 

|2o 

* 

j 

i 

-INTERNAL 

SOURCES 

i 

T777 

y / 

i 

/ / / 

© 

' / / 

c = €0  n2 

♦ 


APERTURE 

PLANE 


E? 

(APERTURE 

FIELD) 


(b) 


T7777  7; 

Figure  la,b 


////  V //////  / 


i m 


Figure  lc,d 


..  ..  U|U  , 1 .. 


9 


n k2 

-J^r  /l(z')M  (z-z')dz'  (1) 

*00 

where  nQ  = ^Jz0  is  the  wave  impedance  of  free  space.  The  kernel 
Mo(z-z')  can  be  given  in  terras  of  its  Fourier  transform  M(ot)  : 

CD 

Mq (z-z ' ) = J Mo(a)exp[ik1a(z-z')]da  (2) 

- 00 

where 

Mo(°0  = J0(AC){?2[H^1)(AC)  - H^tfHOJ^AO]  +J0(AO[P(o)-Q(o)]} 

(3) 


For  thin  wires,  it  is  usually  permissible  to  take  Jq(A£)  = 1,  however, 
if  it  is  desired  to  include  sources  internal  to  the  wire  (such  as  a delta- 
function  generator),  the  product  (Ac) Jq(A£)  must  be  left  intact  in 

order  to  correctly  describe  the  singularity  at  the  source  [14].  The 
first  square  brackets  in  (3)  are  the  contributions  from  the  wire  and  its 
perfect  image;  the  functions  P(a)  and  Q(a)  are  corrections  due  to 
the  finite  conductivity  of  the  earth: 

OO 


P(oO 


_2_ 

iir 


/ 


exp(-2HUj) 

VU2 


dA 


(4) 


Q(ct) 


» exp(-2Hu  ) 

— I 2 ~ dA 

iir  J u_+n  u, 

-00  2 1 


(5) 


In  (3) -(5),  we  have  defined 
ux  = (A2  -c2)i 
= (A2  - A* 


C = (1-aV 


Im(c)  > 0; 


Re(u^)  0 

Re(u2)  > 0 (6) 


s,  ■ (n2-“2)i 


Im(?n)  > 0 
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For  the  integral  (2) , we  choose  the  path  so  that 
C = (1  -aV  (a|  < 1 

= i (a2-l)7  |a|  > 1 

This  is  the  original  contour  shown  in  Fig.  2. 

A second  contribution  to  the  average  E at  the  wire  is  due  to  the 
_ext  -gxt 

external  sources  Jg  and  Jm  (still  in  the  absence  of  the  wire-the  boundary 
condition  at  the  wire  surface  has  not  yet  been  enforced).  Using  the  reciprocity 
principle  [29],  this  "incident"  E (the  first  contribution  may  be  thought 
of  as  "scattered")  can  be  expressed  in  terms  of  the  fields  due  to  a source 
at  the  center  of  the  wire: 

J j j [Eq(x";x)  *JgXt  (*")  - H*(x";x)  • J emXt  (x ')  ] dV  (7) 

V 

-W  — — 

Here  V is  the  volume  containing  the  external  sources  and  Eq(x";x)  and 

H*(x";x)  are  the  fields  produced  at  the  point  x"  due  to  an  equivalent 

(averaged)  z-directed  electric  current  source  Jg  = az<5  (x"-h)6 (y")D (z"-z) 

at  the  center  of  the  wire  x = (h,0,z)  , where  D(z)  is  given  in  Appendix  A 

and  represents  the  effect  of  averaging  a unit  z-directed  current  source 

on  the  wire  surface  around  the  circumference  of  the  wire  [14] 

(so  that  D(z)  -*•  <5(z)  as  A + 0) , E*  and  II*  are  given  for  reference  in 

Appendix  A.  Equation  (6)  allows  us  to  formulate  the  solution  for  arbitrary 

excitation  solely  in  terms  of  the  fields  of  z-directed  electric  currents. 

We  shall  suppose  that  the  boundary  conditions  at  the  surface  of  the 

wire  can  be  put  in  the  form 

K 00 

<EZ>=  2^r  / I(z')Zs(z-z,)dz'  - Veq(z) 


(8) 
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The  term  V6C*  allows  for  the  possibility  of  sources  internal  to  the  wire, 
<Ez>  denotes  the  field  averaged  around  the  circumference  of  the  wire,  and 
the  use  of  the  surface  impedance  convolution  operator  Zs(z-z')  assumes  that 
the  fields  at  the  wire  surface  are  essentially  azimuthal ly  independent 
and  TM  [30].  Hence,  this  boundary  condition  should  be  expected  to  be  valid 
in  many  thin-wire  applications  (see  the  discussion  in  Section  III).  In  this 
manner,  one  can  characterize  a wire  of  finite  conductivity  [10,12,13,19,31] 
as  well  as  a thin  dielectric  coated  wire  (Goubau  line)  [10,19,20,32,33]  or 
a cable  with  a leaky  coaxial  braid  [30,34]. 

Combining  (1)  and  (6)  and  enforcing  boundary  condition  (8)  leads  to 
the  following  integral  euqation  for  I(z): 


oo 

J I (z ' )M(z-z ' )dz ' = -"T -<f  [j  [E;(x";x)  • j^Xt(x")  -H^(x";i)-j^Xt(x")]dV" 


n 

o 1 


♦ V6q(z) 


(9) 


where  x = (h,0,z)  and 


M(z-z')  = M (z-z')  + -j— -Zs(z-z') 

o 1 

This  equation  can  be  solved  by  applying  the  convolution  theorem 


(10) 


f f(z')g(z-z')dz'  = J f (a)g(a)exp(ik1az)da 


(ID 


to  equation  (9)  and  taking  its  Fourier  transform  with  respect  to  z, 
resulting  in 


1(a) 


4 

r^yjfu) 


Jext 

m 


-ik.mz" 

(x')Je  dV" 


(12) 
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where  x^'  is  the  transverse  (to  z)  part  of  x"  = (x^'.z")  and  the  Fourier 


transforms  are  given  by 


1(a)  = — f 1 (z)exp(-ik  az)dz  (13) 

TOO  A 

and  similarly  for  M(ot)  ,Zs(a)  ,Ve<^(a)  ,E^(x^';a)  and  (x’^ ;a)  (the  latter  two  may 

be  found  explicitly  in  Appendix  A)  . The  terrfi  Zs (a) , defined  as  the  Fourier  trans- 
form of  Zg(z),  is  the  axial  impedance  of  the  wire. 

In  the  special  case  when  3ext  = jext  = 0 and  V6C*  = V6  (z) , 

e m 

V n(ot)  = k^/2ir  and  we  recover 

_ or  exp(ik,az)da 

I (z)  = I I (a)exp(ik.az)da  = j (14) 

-°o  ^0  M(a) 

-0Xt  — 0xt 

as  found  in  [24].  Likewise  if  Jg  or  is  set  equal  to  a constant 

vector  times  6(x"-xo),  where  xq  is  some  point  in  space,  we  recover  the 
solution  for  the  arbitrarily  oriented  electric  or  magnetic  Hertzian  dipole 
source  located  at  x . The  expression  for  a VED  agrees  with  the  result  of 
more  direct  derivations  [35,36]. 

Of  particular  interest  to  us  is  the  special  case  mentioned  above 
in  connection  with  Fig.  1(d).  Here 


■=ext  ApO  — p .•  it-,  -?ext  p « 

J = 2E  x a 6 (z  ) ; J = 0 ; V n = 0 
m t z e 


corresponds  to  an  aperture  field  distribution  E^  in  the  plane  z 
For  this  case  (12)  becomes 

i(a)  = nok*MCcO  ff  *z‘  X»o(x,,t:-a)]dS 


where  S is  the  (infinite)  aperture  plane  z = 0.  Inserting  (16)  into  the 


13 


inverse  Fourier  transformation  and  assuming  the  order  of  integration  can 
be  interchanged,  we  obtain 


i(2,  = *kJP‘ 


H*(x";-a)exp(ik  az) 

O L 1 


da  > dS" 


M(a) 


(17) 


Following  [24],  we  extend  the  definitions  of  functions  of  a into  the 

complex  a-plane  as  indicated  following  (5),  specifying  also  0 < arg  W^cir, 

2 9 i 2 i 

where  W = (aD  - a ) , aD  = n/ (n  +1)  because  of  the  presence  of  this 

D D 

third  branch  cut  in  the  a-plane  [17].  For  z>0,  we  deform  the  contour  of 
integration  into  the  upper  half  of  the  a-plane.  The  singularity  structure 

~ - W 

of  the  integrand  is  determined  by  M(a)  (Hq  has  only  branch  cuts  already 

- 2 
possessed  by  M(a)),  which  consists  of  three  branch  cuts  and  a number  of 

3 

poles  in  each  of  the  upper  and  lower  half-planes.  The  contour  deformation 
splits  (17)  into  four  contributions  (if  the  branch  cut  at  a = n is  neglected): 

I (z)  = Ipl(z)  + I 2 (z)  + IB1(z)  + IB2(z)  (18) 


where 


167Ti 


pl » 2 


hi™ 


1 

o 
*— • 

m' 

r 

f 

fa  •< 

'§?  X f 

J z 

\x  ° 

S ' 

i eXpCikianl,2z) 


2 S (19) 


pi  ,2 


H^O^i-a)  n^(xj;-a) 


M1  (a) 


Mi  (a) 

TT 


exp(ikjaz)d^ > dS" 


(20) 


2 

The  contribution  from  the  branch  cut  associated  with  the  square-root 
of  (a2-n2)  is  usually  neglected  under  the  approximation  of  a large  refractive 
index  of  earth,  i.e.  | n | 2 » 1,  although  one  need  not  make  such  an  approximation. 

Usually  only  two  in  previous  investigations  [17].  It  is  seen  in 
section  III  that  the  two  poles  can  degenerate  into  one  double  pole.  This 
special  case  is  handled  in  Appendix  B. 
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Here  a ^ ^ are  the  surface  wave  poles,  the  superscript  1 in  (20) 
indicates  branch  cut  and  the  subscript  tt  denotes  arg  Z,  = tt.  Similarly 
the  superscript  2 in  (21)  indicates  branch  cut  2 and  the  subscript  tt 
denotes  arg  W = n.  The  situation  in  the  complex  a-plane  is  shown  in 
Fig.  2. 


From  the  definition  of  Hq  it  is  clear  that  this  decomposition  is  into 
contributions  from  two  discrete  modal  currents  and  two  continuous  spectra 
of  radiation  currents  (the  branch  cut  integrals.)  Although  the  fields  of 
continuous  spectrum  modes  do  not  satisfy  the  radiation  condition  individually, 
they  are  individually  bounded  at  infinity,  and  when  excited  by  a finite 
source,  they  must  collectively  satisfy  the  radiation  condition.  It  is 
shown  in  Appendix  C that  if  E°  is  replaced  by  the  modal  field  corresponding 
to  a discrete  mode,  then  the  only  nonvanishing  contribution  to  (18)  is  the 
corresponding  discrete  modal  current. 

Now  suppose  E°  were  arbitrarily  close  to  the  E-field  of  one  of  the 
discrete  modes  (which  decays  at  infinity  in  the  transverse  plane).  For 
definiteness,  let  us  say  that  E°  is  the  field  pattern  of  the  discrete  mode 
of  interest,  truncated  at  a circle  of  finite  radius  R,  and  equal  to  zero 
for  p = j x | > R in  the  transverse  aperture  plane  (z  = 0).  E°  then 
differs  from  the  exact  modal  field  by  a quantity  (see  footnote  3)  of  order 
exp(ikj4pR)  or  exp(ikjWpR)  (whichever  is  larger),  where  ^=(1-0^)^  and 
2 2 i 

W^=  (otg  - dp)  define  a transverse  wavenumber  and  a "spreading"  wavenumber 
(describing  the  behavior  of  the  mode  as  y -+  ±°°  along  the  air-earth  interface) 


ORIGINAL 

CONTOUR 


figure  2 

Contour  deformation  in  the  u- plane 
- « ) » lm(W)  > 0;  r.  = (l-a“)a 
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characteristic  of  the  given  mode.  Using  orthogonality,  it  is  then  possible 
to  bound  the  quantities  (20)  and  (21)  by  a value  which  tends  exponentially 
to  zero  (recalling  the  choice  of  Riemann  sheets  for  Wand  C)  as  R + 

Hence  using  only  orthogonality,  without  assuming  completeness  of  the 
modes,  it  has  been  demonstrated  that  either  of  the  discrete  modes  can  in 
principle,  be  excited  with  arbitrarily  high  efficiency  by  taking  a large 
enough  truncated  replica  of  the  modal  E-field  as  the  aperture  excitation. 

It  is  then  apparent  that  given  such  an  excitation,  the  total  current  on 
the  wire  is  essentially  that  of  the  mode  under  consideration  for  wire 
lengths  of  practical  interest. 

III.  Further  numerical  investigation  of  the 
behavior  of  the  modal  propagation  constants 

In  a previous  report  [17],  values  for  the  propagation  constants  of  the 
two  discrete  modes  for  a single  value  of  n and  a and  various  values 
of  h were  presented  for  the  case  of  a bare  wire.  These  had  been  calcu- 
lated numerically  by  using  approximations  of  the  Sommerfeld  integrals 
appearing  in  the  modal  equation  in  terms  of  known  analytic  functions.  The 
approximations  were  based  on  the  assumption  of  a highly  but  finitely  con- 
ducting earth  and  sufficiently  large  height  H.  In  order  to  verify  the 
accuracy  of  the  approximate  modal  equation,  we  compared  the  solutions  of 
the  modal  equation  obtained  from  the  approximate  one,  to  those  obtained 
from  the  exact  one,  wherein  the  Sommerfeld  integrals  were  evaluated  by 
numerical  integration.  More  specifically,  the  "exact"  modal  equation  used 
is  of  the  form  MQ(a)  = 0,  where  Mq  is  given  by  (3).  An  iterative  scheme 
similar  to  the  commonly  used  Newton's  method  algorithm  was  used  to  find 


1 


i 

i. 

I 
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the  roots  of  the  exact  equation.  The  results  for  a single  bare  wire,  to 
seven  digit  accuracy,  are  shown  in  Table  1.  It  is  particularly  gratifying 
to  see  that  the  results  obtained  from  the  approximate  modal  equation  are 
excellent  for  the  largest  h (within  10  , compared  to  those  obtained  from 

the  exact  equation)  and  even  for  the  smallest  h,  still  gives  errors  of 
less  than  10'^.  Within  the  limits  of  validity  outlined  by  these  results, 
it  is  felt  this  comparison  provides  adequate  justification  for  use  of 
either  exact  or  approximate  modal  equations  in  subsequent  calculations. 

It  seems  appropriate  at  this  point  to  discuss  the  only  approximation 
used  in  deriving  the  modal  equation,  i.e.,  that  the  current  on  the  wire 
may  be  taken  to  be  z-directed  and  averaged  over  the  wire  surface,  and  that 
an  average  boundary  condition  around  the  circumference  may  be  used.  This 
would  hold  exactly  for  a symmetrically  excited  wire  in  the  absence  of  the 
earth,  but  due  to  its  influence,  the  so-called  "proximity  effect"  will 
cause  the  wire  current  to  have  higher-order  Fourier  components  in  its 
angular  distribution  around  the  wire,  as  well  as  azimuthally  directed 
components . 

For  a two-wire  transmission  line  with  imperfect  conductors  the  effect 
of  higher-order  Fourier  components  was  taken  into  account  in  a classic 
paper  of  Mie  [37]  (a  more  accessible  source  is  Sommerfeld  [38];  for  a 
generalization  to  arbitrary  noncircular  conductors,  see  [39]).  For  the 
bare,  imperfectly  conducting  wire  over  the  earth,  some  discussion  of  the 
relative  error  in  incurred  by  the  neglect  of  the  first  higher  terms 

of  the  Fourier  series  is  done  in  [12],  but  a more  complete  treatment  has 
recently  been  done  by  Pogorzelski  and  Chang  [40].  For  the  parameter  values 
used  here,  these  errors  can  typically  be  shown  to  be  much  less  than  one 
percent . 
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This  discussion  emphasizes  that  a previous  particular  result  obtained 
by  Coleman  for  the  limiting  case  of  h + 0 [8,15,20]  must  be  understood  as 
a legitimate  limiting  value  only  when  the  ratio  a/2h  is  maintained 
smaller  than,  say,  0.1,  as  the  limit  h 0 is  taken.  It  is  unclear  what 
effect  a small  but  finite  wire  radius  might  have  upon  the  propagation 
constant  of  a wire  located  in  the  interface. 


In  Fig.  3 the  solutions  of  the  bare  wire  modal  equation  MQ(a)  = 0 


are  plotted  for  several  values  of  a and  h,  with  n = 5.3  + i0.95  (the 


loss  tangent  for  this  case  is  about  twice  that  for  the  case  of  Table  I). 


-3  -4 

Between  the  values  a = 10  X and  a = 1.66x10  X , an  interesting 


phenomenon  is  seen  to  have  occurred.  At  some  value  h in  the  neighborhood 
of  0.2X  , there  must  be  a pair  of  curves  for  some  value  of  a between 
those  of  curves  2 and  3 which  actually  have  this  value  of  h in  common. 
This  phenomenon  is  modal  degeneration  [41-43]  which  is  known  to  appear  in 
the  earth- ionosphere  waveguide  [44]  and  in  the  earth-crust  waveguide  [45], 
for  example.  Further  discussion  will  be  reserved  for  the  conclusion. 

Figure  4 presents  data  for  a dielectric-coated  wire  (Goubau  line) 
located  above  the  earth.  Here  the  modal  equation  is  M(a)  = 0,  with  the 


Zs(a)  corresponding  to  (10)  given  by  [10,  19,  20,  32,  33]: 


klno  «r  (i)  Jo(5rB)Hi1)(?rA)‘Jo(5rA)Hi1)ctrB) 

ZS(«)=--T-TH1  — TTi (22) 

n 

r 


J1(5rB)Hi1)CCrA)-J0(CrA)H1C1)(crB) 


where  nr  is  the  refractive  index  of  the  coating  and 


,2  2 i 

£r  = (nr  - a ) 


B = kb 


In  equation  (3)  for  MQ(a).  A must  be  replaced  by  B,  which  is  now 


the  outer  radius  of  the  wire,  while  A is  the  radius  of  the  inner  conductor. 
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figure  3 

Normalized  propagation  constant  a as  a function  of  wire  height  h 
and  earth  refractive  index  n = 5.3  + i0.95. 
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Figure  4 

Normalized  propagation  constant  « of  a dielectric-coated  wire  of 
conductor  radius  a,  coating  radius  b,  coating  index  n , at  a 
height  h above  earth  of  refractive  index  n. 
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Once  again,  we  find  a modal  degeneration,  this  time  occurring  between 
n^  = 1.10  and  1.11  (curves  3 and  4 of  Fig.  4)  at  a height  of  about  0.235A. 

Thus,  Kikuchi's  conclusion  [13,19]  regarding  the  continuous  transition 
from  a TEM-like  wave  to  a surface- like  wave  must  be  modified  somewhat. 

For  the  example  of  Fig.  4,  if  the  coating  index  nr  is  large  enough 
(>  1.11),  then  the  TEM-like  wave  at  low  height  does  continuously  transform 
into  a surface-like  wave  (exactly  a surface-wave  if  h ->■<*>)  with  increasing  h. 
The  second  mode  (an  "earth-attached"  or  "substrate-attached"  mode)  remains 
in  the  vicinity  of  the  branch  point  cig,  only  slightly  influenced  by  the 
proximity  of  the  TEM  mode  at  h ~ 0.24A  . For  smaller  values  of  n^,  on 
the  other  hand,  the  mode  which  is  earth-attached  at  small  h transforms 
continuously  into  the  surface  wave  for  large  h,  while  the  "transmission- 
line" mode  transforms  into  an  earth-attached  mode  for  large  h. 

IV.  Field  plots  for  the  discrete  modes 
A computer  program  to  determine  the  transverse  electric  and  magnetic 

4 

field  distributions  for  each  of  the  discrete  modes  has  been  developed. 

We  present  in  Figs.  5-7  the  relative  magnitude  of  the  electric  field 

for  both  discrete  modes  directly  above  and  below  the  wire  (y  = 0)  as  a 

function  of  the  vertical  coordinate  x.  For  the  cases  n = 1.1,  1.25. 

r 

the  dotted  line  corresponds  to  a mode  which  has  a pronounced  surface  wave 
character,  while  the  "earth  attached"  mode  (solid  line)  decays  more  slowly 
above  the  wire.  For  nr  = 1.0,  neither  of  the  modes  is  distinctly  identi- 
fiable as  surface,  "TEM"  or  earth-attached,  and  the  field  distributions  of 
the  modes  are  somewhat  less  symmetrical. 

— - 

These  are  essentially  those  given  in  Appendix  A. 
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Magnitude  of  electric  field  of  (loubau  line  (in  relative  units)  on 
y = 0 axis.  n=5.3  + i0.45,  h=0.24A,  a = .007A,  b = .01A,  n = 1.0, 
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figure  6 

Magnitude  of  electric  field  of  Goubau  line  (in  relative  units)  on 
y = 0 axis,  n = 5.3  + i0.45,  h = 0.24A,  a = . 007A,  b = .01A,  nr  = 1.1  , 


(o  = 1.0019  + iO. 01132),  (a  = .98755  + iO. 00655b). 
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Figure  7 

Magnitude  of  electric  field  of  Goubau  line  (in  relative  units)  on 
y = 0 axis.  n=5.3  + iO.45,  h=0..24A,  a = . 007A  , b=.01A,  n =1.25, 
(a  = 1.0123  + 0.01134),  __  (a  = .9854  + i0.00577)r. 
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In  the  cross-sectional  plots  of  Figs.  8-13,  the  greater  transverse 
spread  of  the  "earth-attached"  mode  can  be  seen  as  compared  with  the 
transmission- line  mode.  The  fields  are  all  nearly  linearly  polarized 
for  the  latter,  with  the  exception  of  points  very  far  removed  from  the 
wire.  These  plots,  in  conjunction  with  the  excitation  discussion  of 
section  II,  allow  a design  of  an  efficient  experimental  excitation  scheme 
for  selectively  producing  one  of  the  discrete  modes  while  minimizing  the 
other.  This  is  a necessary  development,  if  practical  utilization  of 
these  modes  is  to  be  achieved. 

V.  Conclusion 

In  this  report,  a proof  has  been  given  that  the  current  on  a thin  wire 
parallel  to  a conducting-half  space  may  be  uniquely  and  completely  decomposed 
into  discrete  and  continuous  modal  components.  In  addition,  any  of  the 
discrete  modes  may  be  excited  with  arbitrarily  high  efficiency  (and  to 
the  exclusion  of  the  other  modes)  by  appropriately  matching  an  aperture 
field  to  that  of  the  desired  mode.  This  means  that  even  though  the  con- 
tinuous or  "radiation"  type  currents  must  always  eventually  dominate  the 
discrete  modal  currents  on  a lossy  structure  such  as  this  [28],  the  discrete 
contributions  can  be  made  to  dominate  over  an  arbitrarily  long  portion  of 
the  wire,  so  that  the  concept  of  modal  representation  will  still  have 
important  utility  in  practical  applications.  Moreover,  regardless  of  which 
modal  contributions  (if  any)  are  dominant,  the  modal  representation  is  always 
rigorous,  and  the  attendant  orthogonality  properties  are  useful  in  solving 
problems  wherein  the  wire  is  not  infinite  or  possesses  discontinuities, 
since  such  situations  can  be  treated  by  straightforward  generalizations  of 
closed  waveguide  techniques  [26]. 
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A numerical  investigation  of  the  behavior  of  the  discrete  mode 
propagation  constants,  more  extensive  than  those  which  have  appeared 
previously,  shows  that  the  existence  of  two  discrete  modes  (called  "trans- 
mission-line" and  "earth-attached"  in  [17])  rather  than  only  one,  is  not 
due  to  any  approximation  made  on  the  modal  equation,  but  follows  even 
when  the  exact  integral  representations  of  the  functions  P(a)  and  (ot) 
appearing  therein  are  used.  Moreover,  trajectories  of  the  modal  propagation 
constants  in  the  complex  a-plane  have  been  given,  and  the  possible 
existence  of  modal  degeneration  has  been  demonstrated. 

It  should  be  emphasized  that  degeneracy  of  this  kind  (algebraic  but 
not  geometric,  in  the  language  of  linear  operators  [46])  is  not  basically 
a similar  phenomenon  to  that  which  occurs,  for  instance,  in  closed  circular 
waveguides  between  TE  and  TM  modes.  The  latter  (which  is  both  algebraic 
and  geometric,  and  could  be  designated  a "semi-simple"  mode  following 
Kato  [46])  is  accidental,  a result  of  the  high  symmetry  of  the  waveguide, 
and  in  this  case  there  remain  two  orthogonal  eigenfunctions  (mode  functions) 
corresponding  to  the  same  modal  eigenvalue.  For  the  nongeometric  degeneracy, 
however,  only  one  modal  eigenfunction  remains,  but  the  Green's  function 
(resolvent  kernel)  possesses  a double  pole  at  the  degenerate  eigenvalue. 

Perhaps  the  most  commonplace  example  of  a nongeometric  degeneracy 
(though  one  not  usually  thought  of  from  this  viewpoint)  is  at  cutoff  of  a 
mode  in  a closed  waveguide  [47].  In  this  case,  the  forward-and  backward- 
traveling modes  have  the  same  propagation  constant  (zero)  and  as  a con- 
sequence the  norm 
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can  be  shown  to  vanish,  and  the  excitation  equations  must  be  modified  in 
a manner  similar  to  Appendix  B. 

The  question  of  expanding  the  fields  into  modes  will  be  taken  up 
in  a future  report.  In  particular,  this  will  require  construction  of  the 
dyadic  (tensor)  Green's  functions  as  a generalized  product  of  eigenfunctions, 
as  well  as  a treatment  of  the  generalized  orthogonality  and  completeness 
properties  of  a degenerate  discrete  mode. 
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Appendix  A 


In  this  Appendix  we  present  the  expressions  for  the  fields  Ew  and  HW 

o o 


of  the  z-directed  electric  current  source 
J = a 5(x"  - h)6  (y")D(z"-z) 

6 Z 

where  D(z"  - z)  is  the  generalized  function 

k r°°  ik  a(z"-z) 

D(z"  - z)  = 2 W J J0(AC)e  1 


da 


(A.l) 


(A.  2) 


The  ^0  term  in  (A.  2)  reflects  an  averaging  of  a dipole  source  at  z over 
the  circumference  of  the  wire  (see  [14]).  The  field  expressions  are  found 
in  a straightforward  manner  (see,  for  example,  [15])  and  given  in  terms  of 
the  Fourier  transforms  as  functions  of  x”  = (x",yn): 


_w  _ _ r1*  ~w  _ +ik  a(z"-z)  . 

Eo(x";x)  = j E^(x-';a)e  1 da 


— w — — r - w — i 

Hq(x";x)  = J H“(x^;a)e  1 


+ik  a(z" -z) 


da 


(A.  3) 


The  transforms  are  given  in  terms  of  two  Hertz  potentials  it  and  tt*  : 


3tt 

3tt* 

= ikj 

ax77  + no 

— } 
3y" 

3tt* 

(A. 4) 

Eoy(7t;<0 

= ikj 

<“  IF  • no 

3^} 

Eoz(It:“) 

= (k2 

- k2a2)ir 
1 z 

Hox(*t;<,) 

aft* 

aff 

= ikj 

(a  3F  - (c, 

/V 

1} 
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(A. 5) 

n“zup«) 
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APPENDIX  B 

In  this  appendix  we  give  the  derivation  for  the  excitation  of  the 
discrete  modal  currents  in  the  case  when  the  two  discrete  modes  have 
degenerated  (ctpl  = = ap) . This  analysis  is  along  the  lines  of  those 

in  [41,42]. 

Proceeding  from  (17j,  let  us  find  the  contribution  of  the  double  pole 

at  a . Now,  the  residue  of 
P 


Hw  OL" ; -a)  exp(ik  oa> 
a t 1 


M(a) 


at  a is  given  by 
P 


.•  2 ^ ii  | 

3 j (a-a  ) H*  (xt  ; -a)  expCik^z)) 

— ^ *■  ■ — 

[ M(a)  | a = a_ 


CB.l) 


2 to  [>£  (xt;  -aJexpUk 


M"(a  ) 


la2)] 


2M" 


(°p) 


a=ap  3[m"  (otp)  ] 


Ho  ^xt > 'ap^exP^i^iapz^ 


Thus  the  pole  contribution  to  (17)  is  given  (for  z > 0)  by: 


I (z)  = ~TT 

P n„k 


o 1 M 
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In  addition  to  the  contribution  from  the  ordinary  discrete  modal  field 
H*  (x";-a  ) at  the  aperture,  there  is  a contribution  from  the  field 
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which  does  not  have  the  behavior  of  the  ordinary  mode  (the  fields  will  be 
somewhat  mere  spread  out).  Following  the  nomenclature  of  the  mathematical 
theory  of  linear  operators,  this  may  be  termed  an  "adjoined  mode"  [41],  it 
should  also  be  noted  that  the  excitation  of  the  ordinary  modal  current 
possesses  a factor  growing  linearly  with  z.  This,  of  itself,  does  not 
imply  any  violation  of  conservation  of  energy,  since  a mode  at  degeneracy 
is  self-orthogonal  [41,42]  and  energy  is  readily  exchanged  between  it 
and  the  adjoined  mode. 
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APPENDIX  C 

In  this  Appendix  we  show  that  when  the  aperture  field  is  replaced  by 
the  modal  field  for  one  of  the  discrete  modes,  then  all  the  terms  of  (17) 
vanish  with  the  exception  of  the  corresponding  discrete  modal  current.  These 
arguments  are  based  on  orthogonality  proofs  by  R.G.  Olsen  [48]  which  are 
similar  to  those  given  for  related  problems  in  [49-51]. 

Suppose,  for  definiteness,  that  E°  = , corresponding  to  the  propa- 
gation constant  ot^.  Now  the  field  = H*  (x^.;  -a^)  is  a transpose 

modal  field  (propagating  in  the  -z  direction)  whose  fields  are  related  to 
the  forward  mode  by  (see  A. 4 and  A. 5 in  Appendix  A) 

Et2  ='Et2  ’ Ez2  = +Ez2  ’ “t2  = + E*t2  ’ Hz2  ='HZ2 


Both  discrete  modal  fields  satisfy  Maxwell's  equations  with  no  sources, 
as  well  as  the  boundary  conditions  (in  the  uni form- current  approximation) 
at  the  wire. 

Using  the  Lorentz  reciprocity  relation 


V ■ (E+xll'  - E'  xHp  = 0 

where  V = + ikj(oi  ^ - ap2^az  > we  integrate  over  a cross-section  S 

of  large  radius,  and  apply  the  divergence  theorem  to  the  portion  involving 
V^.  Because  of  the  exponential  decay  of  the  integrand  in  the  transverse 
direction,  the  resulting  line  integral 


j\  ' (q*n‘  - |-xh*: 
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tends  to  zero  as  the  boundary  C of  S recedes  to  infinity.  Thus  in  the 
limit  as  S becomes  the  entire  (infinite)  cross-section,  we  have 

■ <E!XR 2 ' B2*'',PdS  . 0 capl  K °p2)  (c.2) 

In  a similar  manner,  one  shows  that 

* ff2'  • *2  * BPdS  " 0 (“pi  l*-“p2>  < C.W 

combined  with  (A.l)  and  (A. 2),  (A. 3)  shows  that  the  integral  I ^ in  (18) 
vanishes . 

If  it  is  noted  that  all  that  was  used  in  (A.l)  - (A. 3)  was  the  fact  that 
the  fields  satisfied  Maxwell's  equations  and  the  boundary  conditions,  and 
that  the  integrands  decayed  exponentially  at  infinity,  we  can  see  easily 
that  the  integrals  I ^ and  Ig2  must  vanish  in  precisely  the  same  way. 
Thus,  if  the  aperture  field  could  precisely  match  that  of  one  of  the  discrete 
modes,  the  only  contribution  to  the  current  on  the  wire  would  be  that  of 
the  same  discrete  mode. 


